Holonomy invariance, orbital resonances, and kilohertz QPOs 
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Abstract 

Quantized orbital structures are typical for many aspects of classical gravity (Newton's as well as Einstein's) . 
The astronomical phenomenon of orbital resonances is a well-known example. Recently, Rothman, Ellis and 
Murugan (2001) discussed quantized orbital structures in the novel context of a holonomy invariance of 
parallel transport in Schwarzschild geometry. We present here yet another example of quantization of orbits, 
reflecting both orbital resonances and holonomy invariance. This strong-gravity effect may already have been 
directly observed as the puzzling kilohertz quasi-periodic oscillations (QPOs) in the X-ray emission from a 
few accreting galactic black holes and several neutron stars. 

1. Quantization of circular orbits: holonomy invariance 

In a curved spacetime, consider a vector under parallel transport along a closed curve. The initial and final 
direction of the vector are not holonomy invariant: in general they will be different after one loop. But could 
they be holonomy invariant after N loops? Rothman, Ellis and Murugan [1] show that for geodesic motion 
along an r = const circle in Schwarzshild geometry (r is the standard Schwarzschild radial coordinate) the 
holonomy invariance occurs after N loops if and only if 
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where tq is the gravitational radius, N = rim, and 1 < n < m are two integer numbers. It is clear that the 
holonomy invariance condition leads to a "quantization" of orbits. 



2. Quantization of circular orbits: the closure condition 



Consider now a slightly perturbed circular orbit, r + where («r, and t is the standard Schwarzschild 
time coordinate. The equation that governs the perturbation is that of a harmonic oscillator (see e.g. [2]), 

£ + o# = 0, (2) 

where a dot means d/dt, and oj r is the radial epicyclic frequency. For stable oscillations uj 2 > 0, while for 
ui 2 < perturbations are clearly exponential in time and thus unstable. 

With Euclidean geometry of space (which is assumed in Newton's theory), and for a spherically sym- 
metric gravitational potential, orbital frequency and radial epicyclic frequency are always real and equal, 
fl 2 = > 0. Therefore, all slightly perturbed circular orbits close after exactly one loop, and they are all 
stable. This conclusion, of course, holds not only for slightly perturbed orbits which we are now discussing, 
but for all bound, i.e., with negative energy, Newtonian orbits in a spherical potential — they all are stable 
ellipses, as was known already by Kepler. 
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However, ft 2 ^ uj 2 if the gravitational potential is non-spherical, or the geometry is non-Euclidean, which 
means that, in general, perturbed circular orbits do not close after one loop. It may be in addition that for 
a range of radii uj 2 < 0, which means that orbits are unstable in this range. These facts are well-known in 
the context of Schwarzschild geometry, where [2] , 
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and therefore the slightly perturbed circular orbits do not close after one loop (this is the reason for the 
relativistic advance of the perihelion of Mercury) . In addition the perturbations are unstable in the region 
f < r ms = 3r<3, with r ms being known as the marginally stable orbit radius. 

It is perhaps less known, however, that the same effects do happen also according to Newton's theory 
in a gravitational field of a sufficiently non-spherical body. It has been recently demonstrated [3] that the 
marginally stable orbit exists outside rapidly rotating Maclaurin spheroids, as well as outside strange stars 
(quark stars with an equation of state comfortably consistent with the main-stream views on super-dense 
matter) with a very large quadrupole moment, but well in the Newtonian regime (r ms > i?» >> re)- Here 
i?* is the stellar radius. 

Returning to the Schwarzschild geometry, let us state some obvious consequences of Eqs. (3). The 
slightly perturbed circular orbits close after an integral number of loops if and only if nft — muj ri or (which 
is the same) if and only if the closure condition is fulfilled, 
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where 1 < n < m are two integer numbers. 

We see that the closure condition (4) leads to a similar (but different) quantization of circular orbits 
as does the holonomy invariance condition (1). Rothman et al. noticed that holonomy invariance occurs 
at circles outside of the free circular photon orbit, r > r p h = (3/2)rQ. They stressed that this could be 
intuitively understood in terms of the optical geometry, i.e., 3-space conformally rescaled in such a way that 
light trajectories are geodesic in it [4, 5]. The fact, most relevant in the context of holonomy on timelike 
circular geodesies, that r p h coincide with the circle of a minimal circumference length, is an obvious triviality 
in optical geometry. Similarly, the closure condition can only be fulfilled for circles outside the marginally 
stable orbit for test particles, r > r ms = 3ro- The innermost orbit for which the closure is possible has a 
minimal curvature radius, a fact which again is most relevant, and trivial to understand in optical geometry. 



3. Quantization of orbits in accretion disks: resonances and kHz QPOs 



In numerous sources that astronomers directly observe, matter is being accreted onto the central black black 
hole through a vertically thin accretion disk: the vertical thickness H is much smaller than the corresponding 
radius r. For inner parts of accretion disks that are relevant in the present context, typically h — H/r w 0.1. 
In such disks, and in the relevant range of radii (a few times r ms ), matter moves on orbits not much different 
from circular geodesies: pressure, viscosity, radiation, magnetic field, and other effects could be considered 
as being very small perturbations. Therefore, in thin accretion disks, there should be present oscillation 
modes whose frequencies (in the considered range of radii) are very close to cj r , as well as other modes whose 
frequencies are close to the meridional epicyclic frequency ojg, i.e., the frequency of oscillatory motion of 
a test particle about the orbital plane. In Schwarzschild geometry, uog — ft, and in Kerr geometry uog is 
somewhat smaller than ft. The functions ft(r), w r (r), LOQ(r) are explicitly known for both Schwarzschild and 
Kerr geometry (see e.g. [2, 6, 7]). 

The general statement that the frequencies of global modes in thin accretion disks should be determined 
by the space-time metric is well supported by detailed studies of thin accretion disk oscillatory modes, or 
as it is called "diskoseismology" , pioneered by Kato and developed by Wagoner and collaborators [2,6]. It 
has been suggested that prominent g— modes and c-modes can be identified with the "stable" frequencies 
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observed in black hole candidates [8], but the frequencies of these modes are not expected to be in rational 
ratios. Other frequencies have also been investigated [2,6,7,9]. The well-known and important global p-modes 
have frequencies close to w r . Other, local, modes have frequencies close to cue [10]. In these cases, pressure 
and other non-geodesic effects contribute only to a small shift of the two frequencies. 

Our suggestion, is that the two frequencies could be in a rational ratio, ncog — muj r , and when this 
happens, a resonance may occur and the resulting oscillation is prominent (and relatively easily detectable) 
for some, but not all, values of n/m. Two of us have previously pointed out [11,12,13] that such resonant 
frequencies may indeed be observed as kilohertz quasi-periodic oscillations (QPOs) in X-ray emission of 
accreting black holes and neutron stars (for a review of QPOs properties see [14]), and show elsewhere [10] 
that it follows naturally from this rational-ratio constraint and from fundamental properties of the relevant 
nonlinear oscillators that 3:2 and 5:3 ratios are expected among the observed frequencies, as those connected 
with the fastest growing modes. In the black hole source GRO J1655-40 the two high frequencies observed 
are 450 Hz and 300 Hz [15], i.e., they are in a 3:2 ratio [12,13], as expected in our theory. The same 3:2 
ratio of frequencies may have also been observed in the black hole source XTE J1550-564 [16]. In the black 
hole source GRS 1915+105 the two frequencies observed are 70 Hz and 42 Hz [17], i.e., they are in the also 
expected 5:3 ratio [12]. 

We will not further discuss here our theory of excitation of such coherent perturbations in accretion 
disks. We only point out that, most typically, the relative size of the perturbation is of the order of h and 
therefore its relative surface of the order h 2 « 0.01. If the resonant modes are excited, all properties of 
thin accretion disks should oscillate, and in particular there should be an oscillatory part to the emitted 
X-ray flux, of fractional amplitude ~ h 2 i.e., a few percent, as observed. The lifetime of the perturbation 
is regulated not only by shear and dissipation in the disk, but also by the resonance itself, and therefore it 
need not be of the order of the orbital period, but could be much longer. 

Note, that in the case of Schwarzschild geometry one has [2] cog = and therefore the resonances occur 
at the "quantized" orbits of radius given by Eq. (4). For the Kerr black holes loq ^ O. This is good news: 
the Kerr angular momentum parameter a appears directly in theory (governing some of its quantitative 
aspects), and thus could be estimated from fitting theory to observations [13, 18]. 

At present, the QPO phenomenon is the only known strong-field consequence of Einstein's relativity 
that could have been already directly observed (with an adequate precision), and not mixed up with other 
effects that arc difficult to follow in all relevant details. No other observed phenomenon could rival QPOs 
in this respect. In particular, the oft-discussed profiles of the Fe fluorescent line seen in X-ray spectra of 
some active galactic nuclei, and believed to originate close to black holes [19], in fact depend on many "dirty 
astrophysics" factors that are difficult to estimate [20] . The brilliant idea of Narayan, Lasota and others (see 
e.g. [21]) that the observational fact that accreting black holes sources are dimmer than accreting neutron 
stars [22] is a direct evidence for the event horizon, suffers from the same difficulty — unknown factors enter 
relevant details in a way impossible to control today [23] . 

4. Conclusions 

Natural frequencies of oscillation of fluid on nearly circular orbits close to black holes (and neutron stars) 
are a clean property of the strong gravity alone — as is the closely related property of holonomy invariance. If 
Einstein was right, they simply are there, and the theory of gravity need not be supplemented by assumptions 
on the equation of state, turbulence, magnetic field value, opacity, or any other special conditions, in order 
to explain the values of fundamental frequencies in oscillatory behaviour of matter near black holes. There is 
a very serious possibility that the observed QPOs reflect directly these frequencies. The QPO phenomenon 
is unique in its importance for testing strong gravity. 
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